Lower Bound for LMC complexity measure 
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Abstract 

Lower bound for the shape complexity measure of Lopez-Ruiz-Mancini-Calbet (LMC), Clmc, 
is derived. Analytical relations for simple examples of the harmonic oscillator, the hydrogen atom 
and two-electron 'entangled artificial' atom proposed by Moshinsky are derived. Several numerical 
examples of the spherically confined model systems are presented as the test cases. For the homo- 
geneous potential, Clmc is found to be independent of the parameters in the potential which is 
not the case for the non-homogeneous potentials. 
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I. INTRODUCTION 

There are several statistical measures of complexity A given measure becomes 

significant when a rigorous bound on it is known to exist. In this letter, we focus on the 
LMC(Lopez-Ruiz-Mancini-Calbet) complexity [if, Clmc > with the aim of deriving a general 
lower bound for it. The lower bound is tested by presenting (a) the analytical expressions 
for some simple systems: the harmonic oscillator, the hydrogen atom and two-electron atom 
proposed by Moshinsky |3j and (b) the numerical calculations on the spherically confined 
model one and two electron systems {4, 5), 0, Q]. 

Consider a D- dimensional distribution function /(r), with f(r) nonnegative and 

n 

J f(r)dr = 1; r stands for ri,...,r£>. The Shannon entropy [8] and the Shannon entropy 
power are defined as 

S f = - f /(r)ln/(r)dr. (1) 



H f = e s f, (2) 
respectively. The so called disequilibrium D has the form 

D f = J f(v)dr. (3) 

The definition of the LMC complexity measure is 

C LMC = H.D (4) 

r~l 

It is known [9| that the complexity corresponding to probability distributions given by rect- 
angular, triangular, Gaussian and exponential functions in one-dimensional position space 
is given by 1, (2/3)(e 1//2 ), (e 1//2 )/2, and e/2, respectively. The rectangular probability dis- 
tribution, by definition, corresponds to the minimum statistical complexity. We shall now 
derive the lower bound for Clmc corresponding to a given one-electron density. 

II. INEQUALITY FOR THE LMC COMPLEXITY 

To derive a lower bound for the LMC complexity we cite Theorem 2 of the paper of 
Yanez, Angulo and Dehesa [10(]. The position-space entropy S e of an A-electron system in 



2 



a physical state characterized by the (normalized to N) one-electron density g(r) fulfills the 
inequality 

S e +<\ng(r)><NlnU^^, (5) 
where g(r) is an arbitrary positive function. From the relationship between the Shannon 



entropies coming from densities normalized to 1 and N 



10|: 



S f = ^ + \nN (6) 
and taking g = f 2 , we obtain the inequality 

S/ + m(7 / 2 (r)^ >0 (7) 
From the definition of the LMC complexity (jlj) we obtain the upper bound 

\nC LMC >0 (8) 

or 

C LMC > 1. (9) 

In the next section analytical expressions are presented for some simple systems, such as the 
harmonic oscillator, the hydrogen atom and two-electron atom proposed by Moshinsky. 

III. ANALYTICAL EXAMPLES 

Consider first the one dimensional box problem i.e. the infinite square well in one dimen- 
sion: V(x) = 0, if \x\ < a and V(x) = oo, if |x| > a. Even solutions are 

^ n = Acos(^^-J, (10) 



if \x\ < a. Odd solutions are 

ih~ _ = R sin , 

V 2a 



i> ,n = B sin — ), (11) 
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if \x\ < a. n — 1,2, A — B — l/a 1//2 is the normalization constant. The disequilibrium 
takes the form: 

D = 7- ( 12 ) 
4a 

for all eigenfunctions. The Shannon entropy can be written as 

-)du), (13) 



S P = In a — 1 — / In ( cos ^ 



V 2 



for the even solutions and 



S = \na — 1 — J In ^sin ^— ™^jdiij, (14) 
for the odd solutions. Eqs. (JI2"]) . ([TBI and ((111) leed to the expression 

hLC LM c = S + hiD = lnfy-l- J ln(cos(^p)du) (15) 

and 

lnC iM c = S + kiD = ]jiH^ -1- J In (sin (™)^) (16) 

for the even and odd cases, respectively. So the complexity for a particle in a one-dimentional 
box contains, except the quantum number n, no parameter not only in the ground- but 
excited states, as well. It is possible to extend this result for the particle in a spherical box, 
PIASB represented by the radial Schrodinger equation 



d R n i 2 dR n i 



2E- 



dr 2 r dr 
where the radial wave function R n i(r) is given by 



1(1 + 1) 



Rnl = (17) 



R nl = NjtirV^E) (18) 



where ji is the Bessel function of the first kind of order I and N denotes the normalization 
constant. With the Dirichlet boundary condition imposed according to R n i(r c ) = 0, one 
obtains through the condition j/(r c v / 2 _ E) = or (r c y/2E) = Ui t k : the energy levels given by 
Ek,i = ^ttj where u^i denotes the kHh) zero of ji. For the ground state, one obtains for the 
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disequilibrium 



St{2n) - gg _ 0.6720709 



The Shannon entropy can be written as 



S r = ln(87rr c 3 ) - 3 + (20) 

7T 

From Eqs. (JT9l) and (1201) we obtain the relation 

C iM c = 1.3207394, (21) 

which is constant with respect to the radius of confinement r c . A similar analysis in the 
momentum space using the Fourier transform of the position space wave function in the 
case of PIASB leads to Clmc = 1.517215. As our second example, we consider free linear 
harmonic oscillator with potential V = l/2kx 2 . Then the ground-state density has the form 

f? = ^exp(-A; 1 /V). (22) 
We can immediately obtain the Shannon entropy: 

k i "\ i 



S = - l "{^) + 2 < 23 > 
and the disequilibrium: 

k 1/4 

From Eqs. (1231) and ([241 we are led to the relation 

In C LMC = S + ln£> = i(l-]n2). (25) 
Thus the complexity of the linear harmonic oscillator in the ground state is 

C LM c = e^ 1 -^. (26) 
This is a remarkable result as it is a constant, independent of k. 
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The third example is the free hydrogen atom, or hydrogen-like atomic ions. The ground- 
state density takes the form 

Z 3 

g = — exp(-2Zr), (27) 

71 

where Z is the atomic number. The Shannon entropy can be written as 

S = 3-\n (^j, (28) 
while the disequilibrium takes the form: 

Eqs. ( 1281) and ( 1291) lead to the expression 

InCWc = S + lnD = 3(1 -ln2), (30) 

so the complexity of hydrogen-like atomic ions in the ground state is 

C LMC = e^ 2 \ (31) 

Here, again there is no parameter in the result. 

In the Moshinsky model two electrons with antiparallel spins interact harmonically in 
isotropic harmonic confinement. The Hamiltonian has the form 

H = \ (-V? + TD + 1 - + rl) + l -Krl 2l (32) 

where 



ri2 = t x - r 2 (33) 
and K is the coupling constant. Introducing relative (Eq. fl33|) ) and centre of mass coordinates 

R = \ (ri + r 2 ) (34) 
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the Schrodinger equation can be separated leading to the equations 



2 

and 



- V 2 , + 2R 2 ) m UM = E VM m VM (35) 



\ ("V r 2 12 + \{\ + K)r\ 2 ) ^ RM = E RM ^ RM . (36) 



The ground-state density of the Moshinsky atom has the form llj 



g(r) = aexp (-br 2 ), (37) 



where 



^ \ 3/2 

a = 2 ( -J , (38) 



2a - 1 

b = (39) 



and 



« = ^(1 + (1 + 2K) 1/2 ). (40) 

From Eqs. (JTJ) , (J3]) , (SJ) and ( 1371) we arrive at the expression 

S + \nD= -(l-ln2). (41) 
2 

Consequently, the complexity of the Moshinsky atom in the ground state is 

C LMC = e^ 2 \ (42) 



This is again remarkable as it is a constant, independent of the coupling constant. In 
all the analytical results treated in this section, the parameter-free result for Clmc is is 
consequence of the homogeneous property of the corresponding potential. This has been 
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proven earlier in the cases of the composite information measures 12j and for complexity 
6|. When the homogeneous property of the potential is lost, for example, under the 

confined conditions as we shall see in the next section, the Clmc values become dependent 

on some parameter of the potential. 



IV. NUMERICAL RESULTS ON CONFINED MODEL SYSTEMS 

The spherically confined Coulomb potential was originally introduced 4j in order to study 
the effects of high pressure on the dipole polarizability of hydrogen atom. In this model, the 
confining potential is defined in terms of an impenetrable spherical cavity of a certain radius 
of confinement, R. The spherically confined Coulomb potential for hydrogen-like atoms, for 
example, is given by 

V(r) = for r < R, 

r 

= 00 for r > R. (43) 

We have carried out the calculations of Clmc corresponding to the following spherically 
confined model systems (a) particle in a spherical box (PIASB) (b) spherically confined 
isotropic harmonic oscillator (SCIHO),(c)spherically confined hydrogen atom (SCHA), all 
in position and momentum space and the two electron systems of (d) a pair of electrons and 
(e) He atom in position space. For the two electron systems we have employed the accurate 
23-term and 70-term Hylleraas wave functions for the (d) and (e) cases, respectively. The 
numerical procedure adopted by us has been described in the literature [5|, |6|, |7| and will not 
be repeated here. Our results on the calculations of Clmc f° r systems (a)-(d) as a function 
of the radius of confinement, R, are displayed in Fig. 1. 

For the PIASB cases (position and momentum space) , the homogeneous character of 

n n 

the potential leads to a constant value of Clmc [6|, 1121]. In all other cases, the absence of 
homogeneous property of the potential gives rise to a nonlinear dependence of Clmc with 
R. The range of confinements considered by us presents a good test of the lower bound 
on Clmc obtained in Eq.(9). As is evident from Fig. 1, this bound is obeyed for all R 
values across a diverse set of quantum confined systems. The general lower bound on Clmc 
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derived in the present work is a significant addition to the characteristics of this statistical 
measure. 



V. CONCLUSIONS 



In summary, we derived a general lower bound for the shape complexity measure of 
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9J We also calculated the LMC measure for simple model systems: the 
harmonic oscillator, the hydrogen atom and two-electron 'entangled artificial' atom proposed 
by Moshinsky. We obtained the remarkable result that the LMC measure is constant, that 
is, it does not depend on any parameter. This is characteristic of the homogeneous property 
of the corresponding potentials. A representative set of confined Coulomb potentials have 
been used to numerically test the general validity of the lower bound obtained on Clmc- Due 
to the non-homogeneous property of the potentials in these cases, Clmc becomes dependent 
on the parameters of the potential. 
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FIG. 1: Variation of Clmc corresponding to the spherically confined model systems in position 
and momentum space for (a) particle in a spherical box (PIAB) (b) isotropic harmonic oscillator 
(SCIH0),(c) hydrogen atom (SCHA), and in position space for (d) a pair of electrons and (e) He 
atom. For (d) 23-term , and for (e) 70-term, Hylleraas wave functions have been employed. 
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